Probing spin and orbital Kondo effects with a mesoscopic interferometer 
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We investigate theoretically the transport properties of a closed Aharonov-Bohm interferometer 
containing two quantum dots in the strong coupling regime. We find two distinct physical scenarios 
depending on the strength of the interdot Coulomb interaction. When the interdot Coulomb inter- 
action is negligible only spin fluctuations are important and each dot develops a Kondo resonance 
at the Fermi level independently of the applied magnetic flux. The transport is characterized by 
the interference of these two independent Kondo resonances. On the contrary, for large interdot 
interaction, only one electron can be accommodated onto the double dot system. In this situation, 
not only the spin can fluctuate but also the orbital degree of freedom (the pseudo-spin). As a result, 
we find different ground states depending on the value of the applied flux. When tf) = n (mod 2n) 
{(j) = 27r"l>/$o, where $ is applied flux, and "l>o = h/e the flux quantum) the electronic transport can 
take place via simultaneous correlations in the spin and pseudo-spin sectors, leading to the highly 
symmetric SU(4) Kondo state. Nevertheless, we flnd situations with (p > (mod 2tt) where the 
pseudo-spin quantum number is not conserved during tunneling events, giving rise to the common 
SU(2) Kondo state with an enhanced Kondo temperature. We investigate the crossover between 
both ground states and discuss possible experimental signatures of this physics as a function of the 
applied magnetic flux. 



PACS numbers: 72.15.Qm, 73.63.Kv, 73.23.-b 



I. INTRODUCTION 

Progressive advance in nanofabrication technology has 
achieved the realization of tiny droplets of electrons 
termed quantum dots (QD's) with a high-precision tun- 
ability of the transport parameters.^ One of the most ex- 
citing features of a QD is its ability to behave as a quan- 
tum impurity with spin l/2mi^ At temperatures lower 
than the Kondo temperature (Tk), the localized spin be- 
comes strongly correlated with the conduction electrons 
and consequently is screenedi ^i'^i'^ As a result of the in- 
creasing rate of scattering there arises a resonance at the 
Fermi energy (Ep) in the density of states (DOS) of the 
QD.^ The transmission through the quantum dot is then 
almost perfect. This is the so-called unitary limit where 
conductance reaches 2e^ /h. ^iLi Among many of the ad- 
vantages offered by QD-based devices, we highlight the 
possibility of studying the Kondo effect out of equilibrium 
by applying a dc bias^i or a time-dependent potentialiSiiS 

A natural step forward is the understanding of the 
magnetic interactions of two artificial Kondo impuri- 
ties,"i^^i"i^^i^^i^^i^^i^^i^^. The investigation of double 
QD's is mainly motivated by the possibility of their ap- 
plication as solid-state quantum bits, by using either 
spin or charge degrees of freedomiiSiiSiSfliSi When the 
two QD's are interacting, the orbital degrees of free- 



dom come into pl ay as a pseudo-spin, as shown exper- 
imentally in Ref. '22 23', which may give rise to exotic 
physical scenarios. ^''■^''^ Thus, in a double QD it is pos- 
sible to tune appropriately the gate voltages in order 
to find two charge states almost degenerate. If the in- 
terdot Coulomb interaction is large enough these two 
states are {ni = 1, n2 = 0} and {ni — 0, n2 ^ 1}) 
where fii(2) = ('^i(2)) is the charge state in the dot 
"1" ("2"). This is one of the basic ingredients to ob- 
serve Kondo physics: the existence of degeneracy be- 
tween two quantum states. The Kondo effect is then 
developed to its fullest extent in the pseudo-spin (or- 
bital) sector. We define the pseudo-spin T as follows: 
it points along +(— )z when the electron is at the "1" 
("2") dot: f ^ = (1/2) (ni - na). The pseudo-spin of the 
double QD system can be 1/2 or —1/2, which is quenched 
(screened) via higher order tunneling processes produc- 
ing the so-called orbital Kondo effectw^^^i^ Other real- 
izations of such exotic "orbital Kondo effects" have been 
recently proposed in different QD-related structures as 
.^gjj|25j26j27 w^jigQ ^}jg intradot Coulomb energy for each 
dot is large, then each QD also behaves as an magnetic 
impurity and the conventional Kondo effect is also ob- 
served in the spin sector (5^ = ±1/2). The quantum fluc- 
tuations between these four states [S^ — ±1/2 = {t,i} 
and = ±1/2 = {tT,!!}] lead to an unusual strongly 
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correlated Fermi liquid state in which the (real)-spin and 
pseudo-spin are totally entangled. — In contrast to com- 
mon spin Kondo physics observed in QD's, this new state 
possesses a higher symmetry, SU(4), corresponding to 
the total internal degrees of freedom of the double QD: 
{TtTj itfj i^}- The screened magnitude is now the hy- 
perspin M = Y^aA^" + 1/2) (T*" -1-1/2). Importantly, 
the associated Kondo temperature y^^^'*) jg much higher 
than in the common spin-1/2 Kondo effect in a QD, 
which makes the observation of this spin and pseudo- 
spin entangled state more accessiblei^ Strong entangle- 
ment of charge and spin flip events is also possible in a 
single-electron box (metallic grain) coupled to a lead via 
a smaller quantum dot in the Kondo regime^SiSi Here, 
the spin Kondo physics stems from the screening of the 
spin of the small dot while the pseudo-spin Kondo physics 
emerges when charging states of the grain with (charge) 
Q — and Q = e are almost degenerate. 

The most prominent feature of the Kondo effect is the 
phase coherence experienced by the electrons that par- 
ticipate in the many-body correlated state. Therefore, it 
is thus of great interest to have access to the phase of the 
transmission amplitude in order to give a fully charac- 
terization of the transport properties. The widely known 
Aharonov-Bohm (AB) effect^^ provides us a valuable tool 
to investigate quantum coherence of electrons. When the 
coherence of a circulating electron wave packet enclosing 
a magnetic flux $ is preserved, the result is an extra flux- 
dependent phase shift (0). In the simplest realization of 
an AB interferometer, an incoming electronic wavefunc- 
tion splits into two paths, which join again into the out- 
going electronic wavefunction. Applying a magnetic flux 
which threads this closed geometry, the outgoing wave- 
function acquires a flux-dependent phase, (j) = 27r$/$o, 
where ^ = B/S is the flux, B is the applied magnetic 
field, S is the enclosing surface, and $0 = h/e the flux 
quantum. As a consequence, the transmission is a peri- 
odic function of (f). 

In this work, we consider a double quantum dot embed- 
ded in a prototypical mesoscopic interferometer thread- 
ing a magnetic flux see Fig. Our motivation to 
investigate this system is twofold: (i) there are striking ef- 
fects, such as Fano resonances, which arise already in the 
noninteracting case^MiiS^ and, more interestingly, (ii) as 
the interdot interaction gets stronger, the local density 
of states on the double QD changes drastically.^'^ Here, 
we provide a unified picture of the combined influence of 
wave interference, Kondo effect, and interdot interaction 
on the electronic transport through a double QD in and 
out of equilibrium. 

As we anticipated, the physical scenario in our setup 
will depend much on the strength of the capacitive in- 
terdot coupling between the two dots. When the inter- 
dot Coulomb energy is negligible each QD can accommo- 
date one electron and both spins become screened. We 
find that each QD develops a Kondo resonance at the 
Fermi level Ep = 0. Their interference causes a very 
narrow dip in the differential conductance = dl/dV^c 




FIG. 1: Sketch of the Aharonov-Bohm interferometer con- 
taining two quantum dots attached to two leads. The arrowed 
straight line indicates the tunnel coupling and the dashed line 
represents the interdot Coulomb interaction. 

except at w (mod 27r). In the limit of strong inter- 
dot Coulomb interaction there are two degenerate charge 
states described by the pseudospin = ±1/2. The 
presence of flux allow us to explore two interesting situa- 
tions, (i) when the pseudo-spin is a good quantum num- 
ber [(p = TT (mod 27r)] and (ii) when the pseudo-spin is 
not conserved during tunneling. In the former case the 
SU(4) Kondo state is fully developed whereas far away 
from this symmetry point the conventional SU(2) Kondo 
physics arises. In addition, we will show that shows 
a zero bias anomaly (ZBA) instead of a dip when the 
interdot Coulomb energy is large, which is suppressed as 
(j) enhances and eventually disappears at the destructive 
interference condition [when cf) ~ n (mod 27r)] resulting a 
complete suppression of the tunneling current. Neverthe- 
less, this fact does not prevent to us to observe the highly 
symmetric SU(4) Kondo state since it survives even away 
from (f> — TT (mod 2tt) where the differential conductance 
is not totally suppressed. 

This work is organized as follows: we begin in Sec. ^ 
presenting the theory to treat both limits for the inter- 
dot Coulomb interaction using different theoretical tech- 
niques. We derive the transport properties as well. In 
Sec, mil we present our numerical results and their inter- 
pretation. Finally, we shall end up by summarizing our 
main conclusions in Sec. II VI 

II. THEORY 

The system that we consider is depicted in Fig. ^ 
It is a closed-geometry AB interferometer, where elec- 
trons emitted from the leads are never lost in surround- 
ing gates. Electrons traveling through the device have to 
go either through the upper dot or through the lower dot 
before being transmitted into either the left or the right 
electrode. The enclosing area by the two paths is pene- 



3 



trated by a flux The two reservoirs are Fermi seas of 
electrons described by the Hamiltonian 



■ytn — 



(2.1) 



l=L,R. kM 



where c|^(^) (cz,(fl),A;.(T) is the creation (annihilation) 
operator for an electron in the state k with spin a in the 
lead L{R). The isolated dots are described by J^: 



1=1,2 



Ui2nin2 . 



(2.2) 

The operator dl ^ {di,a) is the creation (annihilation) 
operator, Si is the level position, Ui are the intradot 
Coulomb interaction, and rii.o- — d\ ^di^a is the occu- 
pation number on the dot i. U12 denotes the interdot 
Coulomb interaction between the dots. The tunneling 
between the dots and the leads is modeled by J^t'- 



E 

3 = 1,2' 



■ L,R k,cr 



h.c . 



(2.3) 



The tunneling amplitude Wej in Eq. (|2.3|l from the dot j 
to the lead i is modulated by the external flux <I> thread- 
ing the loop (Fig. ^| and given by 



Wl.2 = 
Wr^2 



(2.4) 



where Vej is the amplitude in the absence of the flux and 
(j) = 27r$/$o with $0 being the flux quantum ($0 = h/e). 
Then the total Hamiltonian is .^^otai = J% + -^d + -^t- 
To make the physical interpretations of our results 
more clear, we perform a few simplifications. First of 
all, we assume identical dots and symmetric junctions; 
i.e., £1 = £2 = ed, Vi = U2 = U, and 14, 1 = Vl,2 = 
Vr,! = Vr,2 = V. This is only for the sake of simplicityiM 
Furthermore, we consider the wide-band limit, in which 
the couplings are independent of energy. Then, the hy- 
bridization of the dot levels with the conduction band is 
well characterized by the parameters 



or, in the matrix notation. 



rL(0) = rz 



1 e 



i<h/2' 



Tr = , 



(2.5) 



(2.6) 



where F^ = irpiV^ with pi being the DOS in the lead £ 
at the Fermi energy (p^ = pR — po). 

Since we are interested in Kondo correlations^, we 
shall mainly concentrate on the Kondo regime [intradot 
charging energy U ^ 00 and localized level —£d ^ 
(Fi + Tr)] for which the fluctuations of the charges in 
the single dots are highly suppressed. For the interdot 
Coulomb interaction U12, we will investigate two opposite 
limits, namely, (i) U12 = and (ii) U12 = 00. In the for- 
mer case, each dot is singly occupied ((^i) = (712) ~ 1) 



and behaves as separate magnetic (Kondo) impurities. 
In the latter case, the double quantum dot system con- 
tains just one electron ((ni -I- 712) = 1). These two limits 
induce striking differences between the resulting Kondo 
effects. Moreover, the interference modulated by the ex- 
ternal flux (p threading the AB geometry leads to an even 
richer variation of the Kondo effects in either case. Our 
goal is to investigate thoroughly these scenarios. For this 
purpose we employ different techniques: scaling analysis 
(valid for T ^ T/f), the slave-boson mean- field theory 
(SBMFT, for T <C TV), and the numerical renormaliza- 
tion group (NRG) method. We elaborate below on these 
approaches. 



A. Scaling analysis 

We derive effective Hamiltonians in the Kondo regime 
for the two limiting cases (U12 = and U12 = 00) and 
discuss their qualitative features at equilibrium by means 
of the scaling theory. 



1. Case U12 

Firstly, we discuss the large capacitance limit between 
the two dots {U12 — *■ 0). As we mentioned, when U12 is 
vanishingly small (and yet Ui,U2 00), the two dots 
are both singly occupied: (ni) = (712) ~ 1 and each 
dot can thus be regarded as a magnetic impurity with 
spin 1/2. In this situation we notice that it is convenient 
and provides a more transparent picture of the system to 
perform the following canonical transformation 



Cl,k,cr 


1 


g+-i7r/4 


C2,k,cr 


"71 


g-i7r/4 



CL,k,a 
CR,k,<7 



(2.7) 



Under this transformation, the Hamiltonian for the leads 
Eq. H2.1|l is rewritten as follows 



^0 = 51 E^'^^Lfc,'^^^^'^''^' 

H=l,2 k,<T 

while the tunneling Hamiltonian Eq. H2.3|l reads 

-^T= EE ^f^-i'^l.k.ad^a + h.C. , 

i=l,2 Ai=l,2 fe,(T 



(2., 



(2.9) 



where 



1^1,1 = V"2,2 = cos ■ 
Vi,2 = V2S = cos ■ 



• — TT 



(2.10) 



Now, the Schrieffer- Wolff transformationi^ of the Hamil- 
tonians Eqs. (|2.8|) . H2.9|) . and H2.2|l leads to the Kondo- 
like Hamiltonian 
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=^^ondo = =^ + -^Jl (Si 



V4(0)<tVi(0)+^|(0)(tV^2(0) 



+ 4^2 (Si 



V;I(0)(tV2(0) + ^|(0)^V^i(0) 



+ ^ J3 (Si - S2) • [^t(0)<T7Ai(0) - V'1(0)*tV2(0) 



-J4Si-S2. (2.11) 



r 



In Eq. H2.11|) , we have adopted the spinor representations 



(2.12) 



(j = 1, 2 and ^ — 1,2) according to which the spin oper- 
ator on the dot j is given by 



(2.13) 



where cr denotes the three Pauh matrices. 

The couphng constants Ji (« = 1, • • • , 4) in Eq. (12.11(1 
are given initially (in the RG sense) by 



Ji = 2^ 



J2 = Ji cos((/)/2) , 

J3 = Ji sin((/./2) , (2.14) 



where jV is the spin degeneracy. Under the rcnormal- 
ization group transformation^'', these coupling constants 
scale as 



^-Po(Ji 



^-2poJiJ2, 

^ = 2poJiJ3, (2.15) 



J 



where f = — log D indicates the renormalization steps {D 
is the bandwidth). J4 is given 



Ji ~ 2poJi2(0)T(D)(l + cos( 



(2.16) 



where T{D) is order 1. J4 corresponds to a ferromagnetic 
RKKY coupling between the spins in the dots. 

Under the renormalization group transformation all 
the system flows to the strong coupling fixed point with 
the ratios J2/J1, J3/J1, and J4/J1 remaining constant. 
In particular, the solution for the initial conditions ((2.14(1 
satisfies the simple properties 



^ = cos(<^/2) , ^ = sin(^/2) , 
with Ji — !• cx) according to the equation 



(2.17) 



^=2poJi^ (2.18) 
From Eq. ((2.17(1 one can easily see that system behaves in 
distinctive ways for different values of flux cf), especially, 
for = (mod 27r) and (f) — ir (mod 27r). In the absence 
of the external flux (0 = 0), J3 = while Ji = J2 = J 
and J4 = I = pqJi/2. Thus the Kondo-like Hamiltonian 
((2.11() is reduced to 



=^Kondo - ^0 + ^ J (Si + S2) • [V'i(O) + 1^2(0)]^ a- [^i(O) + 1/^2(0)] - /Si • S2. 

I 



(2.19) 



This is the two-impurity (characterized by the two spins 
I Si and IS2) Kondo model coupled to a single con- 
duction band (characterized by ipi + ip2 or equivalently 
ci.k,a + C2,fc.cr)- The two spins are coupled to each other 
ferromagnetically (-1 Si • S2, with / > 0). Due to the 
ferromagnetic coupling and to the fact that both spins 
are coupled to the same conduction band, the total spin 
is underscreened at T ^ OjSS Note that a strong RKKY 
interaction may arise from our peculiar geometry since 
both QD's are directly connected to a single channel in 
the leads. Nevertheless, in an actual experimental situ- 



ationSi^S the QD's are far apart and the RKKY inter- 
action may be negligible. Furthermore, slightly above 
= (mod 27r), even for a large ferromagnetic coupling 
|/| 3> Tk = Dexp {—l/2poJ), the spins of the dots added 
in a S' = 1 state become effectively screenedi^ 



For the flux 4 
J2 = while Ji 



= vr (mod 2tt), the coupling constant 
J3 = J4/2 = J. Then the Kondo-hke 
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Hamiltonian H2.11|l is reduced to 

^Kondo = + ^JSl • Vl(0)<T^l(0) 

+ i JS2 • V^(0)crV'2(0) ~ i JSi • S2 . (2.20) 

This model is clearly distinguished from the one in the 
previous case of = (mod 2-k) cf. Eq. H2.19|l . The 
two impurity spins, |Si and -182, of magnitude 1/2 are 
coupled to two independent conduction bands, i/'i and 
'!/'2 (or equivalently ci^fc.(j and C2,fc,cr), individually. The 
ferromagnetic coupling in Eq. H2.20|l does not play any 
significant role in this case, because its coupling strength 
(/) is the same as the exchange coupling between the 
localized spins and the itinerant spins. Therefore, the 
model (|2.20|l corresponds to the usual single-channel spin 
1/2 Kondo model. 

The coupling constants scales according to the renor- 
malization group equation 

^-2poJ^ (2.21) 
and the Kondo temperature is given by 

T.^i.exp(-^). (2.22) 

In the general case ((/) 7^ 0,7r), the two localized spins 
■|Si and ^82 are coupled to two conduction bands t/'i 
and '!/'2, let alone the ferromagnetic coupling with each 
other. Unlike the previous, special case of = tt, the two 
conduction bands are not independent any longer; see 
Eq. H2.11|l . This fact makes the physical interpretation of 
the model rather involved. However, the renormalization 
group flow [see Eqs. p. 171) and H2.18|l ] and the results 
from the numerical renormalization group method (see 
below) suggest that /or any finite flux ((f) ^ 0), the two 
localized spins are fully screened out at zero temperature. 

2. Case C/12 00 

We now investigate the limit of U12 —> 00 where the 
system properties change completely. In this case, only 



one electron is acommodated in the whole double QD 
system, i.e., {rii + 77,2) ~ 1 having either spin t or spin 
|. The orbital degrees of freedom (pseudo-spin) play as 
significant a role as the spin, and the double QD behaves 
as an impurity with four degenerate levels with differ- 
ent tunneling amplitudes depending on the applied flux. 
Due to the orbital degrees of freedom involved in the in- 
terference, the symmetry of the wavefunction is crucial. 
Therefore, in this limit, it is more useful to work with 
a representation in terms of the symmetric (even) and 
antisymmetric (odd) combinations of the localized and 
delocalized orbital channels lii 

In accordance with these observations, we take the fol- 
lowing canonical transformations: 





1 


'1 1 ■ 




dl,cr 


id 


~ 71 


1 -1 




_d2,a_ 



(2.23) 



for the QD electrons, and 





1 


'1 1 ■ 








"71 


1 -1 







(2.24) 



for the conduction electrons. 

Then we identify the pseudo-spin up (down) as the 
electron occupying the even (odd) orbital. More explic- 
itly, taking the four-spinor representation 

V'd = [de,h^ed^d,o^T,do,i] , (2.25) 
the spin and orbital pseudo-spin operators are given by 

^S = ^4rr^,, ^T^Jv^t^Vd, (2.26) 

respectively, where cr (r) are Pauli matrices operating 
on the spin (pseudo-spin) space. Notice that in this 
even/odd basis the dot pseudo-spin has been rotated: 
^ j^!/ ^ _rpy _^ J.X whereas the spin 

remains invariant. 

In terms of the new operators de,^, do.cr, Ce^k,a- and 
Co_fc_cr, the total Hamiltonian, Jiftotai is rewritten as fol- 
lows 



J2 ^fc.'^^".^.'^ + 



Q— e,o k,<T 



''a, I 



+ -{U + 3Ui2)neno 



C/12) {dlade) ■ {diado) -\{U^ C/12) idl^dl^do,idoA + h.c) + Y.^<- ''l,k,sd<.,s + h.c. , (2.27) 



a— e,o k,s 



where Therefore, the even (odd) orbitals are coupled only to 

Ve = 2V cos(0/4) , Vo = 2V sin(0/4) . (2.28) 
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the even(odd)-symmetric combinations of the conduction 
bands. 

To examine the low-energy properties of the system, we 



J 



obtain for all values of 4> the following effective Hamilto- 
nian by performing a Schricffer-Wolf transformation^: 



(2.29) 



r 



where is the first term in Eq. H2.27(l and ip — 

[ipei , ipei , ipo]- , ''Pol] is the spinor of the itinerant electrons. 
Here the effective coupling constants Ji {i — 1, . . . ,Q) are 
initially (in the RG sense) given by 



Jj^^ Js^ 2.yV 



\vl 



J2 = Ji sm{(f>/2) , 

J4 = Ji cos(0/2) , (2.30) 



and scale according to the RG equations (up to the sec- 
ond order in J's) 

^ = 2po {J! + Jl + Jl) , ^ = po J2(3Ji + J3), 



J5 is given and hy^^t^^ 
J5 = 4po|l^|'cos(0/2)ln 



(2.31) 



Sd + D 

Sd-D 



+ S[poJi{l = Q)fT{D) cos(0/2) . (2.32) 

As one can see from the RG equations (|2.31ll . in gen- 
eral, each coupling constant in Eq. H2.29|l scales differ- 
ently under the renormalization group transformation for 
typical behaviors of the solutions at different values of 
flux (p. Importantly, we show now that the system ex- 
hibits a crossover from 0-flux to tt-Aux. Near the 0-flux 
[(p ~ Q (mod 27r)], the double QD odd orbital is com- 
pletely decoupled from the odd-symmetric lead and only 
the even orbital is coupled to the even-symmetric conduc- 
tion lead with K = 2V [see Ea. l|T^ ]. Equation 
then reduces to a model involving only the spin in the 
even orbital |Se (not |S) 



(2.33) 



where J — 2\V\'^/\ed\- This model was already analyzed 
in Ref. l4ll where it was shown that the ground state 
corresponds to a Fermi liquid state with a greatly en- 



hanced Kondo temperature 



,SU(2) 



(due to a coupling doubling of the even orbital to the 
even-symmetric conduction lead, i.e., Vg = 2V) and the 
orbital pseudo-spin gets frozen completely, T =1. (J5 
does not flow to the strong coupling regime). One can 
easily see that the model (I2.33|) is equivalent to the 
twofold orbitally degenerate Anderson model described 
by the common SU(2) Kondo physicsM' 

Near the 7r-flux [0 w tt (mod 2tt)] the exchange cou- 
plings are J = Ji = J2 = J3, and J4 = J5 = 0. The 
corresponding Kondo-like Hamiltonian reads 



j^j^ = i[S. (V/cn/;) + (V'VV) • T 

+ S • (V'VrV') • T] . (2.34) 

This is the celebrated SU(4) Kondo model, where the spin 
and the orbital degrees of freedom become entangled due 
to the third term in Eq. (|2.34() . The RG equation reads 



dJ/d£ = ApoJ^ , 



leading to 



t; 



.SU(4) 



K 



i?exp(-l/4poJ) 



(2.35) 



(2.36) 



Dexp(-l/4poJ) 



As the flux departs from tt, the degeneracy of the even 
and odd orbitals is lifted and the SU(4) symmetry is bro- 
ken, much like a single Kondo impurity in the presence 
of a Zeeman splitting!^ The crossover from the SU(4) 
to the SU(2) Kondo model occurs at a given critical flux 
(pc- From our NRG calculation (see below) we estimate 
(pc « 0.757r. 

This discussion demonstrates the existence of high- 
symmetry Kondo states in double quantum systems with 
interdot interaction in the presence of an Aharonov- 
Bohm flux. We have shown that the magnetic flux criti- 
cally alters the properties of the ground state, resulting in 
a smooth transition from SU(2) to SU(4) Kondo physics. 
Below, we prove that the differential conductance would 
indicate the principal features of this effect. This is im- 
portant since it would serve as a means of experimental 
detection. 
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B. Slave-boson mean- field theory 

In this section, we adopt the so-called slave-boson 
mean-field theory which captures the main physics of 
the Kondo problem^ at sufficiently low temperatures 
(T < Tk). The SBMFT corresponds to the leading 
order in a ^-large expansion, where ^ is the degen- 
eracy of each site. Such a SBMFT has been recently 
applied to study the Kondo effect in noncquilibrium sit- 
uationsi^iiii^ and in double QD's systemsiSJ^. 



impurity Anderson model (J^otai) in terms of the slave 
boson operators. This way the fermionic operator of 
each dot is written as a combination of a pseudofermion 
and a boson operator: di^^- = bifida, where fi^a is the 
pseudofermion which annihilates one "occupied state" in 
the ith dot and b] is a boson operator which creates an 
"empty state" in the ith dot. We include two constraints 
to prevent double occupation in each QD in the limit 
U 1,11-2 — + oo by using two Lagrange multipliers Ai,A2. 
Thus, the Hamiltonian in the slave boson language reads: 



1. Case Ui2 



Firstly we consider the case of vanishing interdot 
Coulomb interaction Uu = 0. We express the two- 



1=1,2 o- 1=1,2 e=L,R k,a 

+ E ^' (E4-A- + ^I^'-i) ' (2.37) 

1=1,2 \ a I 
I 



where W = Wl,^^^J^ . The hallmark of the SBMFT 
consists of replacing the boson operator by its classical 
(nonfluctuating) average. hi(t)l^rJ^ {bi)/V^ = bi 
thereby neglecting charge fluctuations in each dot. This 
approximation is exact in the limit .yV oo, and it cor- 
responds to ^(1) in a 1/yy expansion. At zero tem- 
perature T — 0, it correctly describes spin fluctuations 
(Kondo regime). Then, the mean field Hamiltonian is 
given by 

J^F = + E E ^hi^flafho- 

+EE + H 

i i,k,a- 

+ ^A, , (2.38) 

i 

where Wi^i = biW We obtain a quadratic Hamiltonian 
containing four parameters {bi^2 and renormalized levels 
£1,2 = Ed + Ai,2) to be determined from mean field equa- 
tions These mean-field equations are the constraints 
for the dot z = 1, 2: 

Y.^fl{t)hAt))+-^\h? = ^^ (2.39) 

(T 

and the equations of motion (EOM) of the boson fields: 

E ^^-^ {clA^)hAt)) + = . (2.40) 

e,k,<7 



The next step is to write these mean field equations 
in terms of noncquilibrium Green functions. The lesser 
dot-dot Green function is (z G 1,2) Gf^(t - t') = 

^i{fi<y{'^')fi,'y{^))^ and the corresponding lesser lead-dot 
Green function is G<„.j,j^ ,^{t - t') = -i(c[fc,^(i')/»,^(*)>- 
By applying the equation-of-motion (EOM) technique 
and the analytical continuation rulcs^^''*^ we can relate 
the lesser lead-dot Green function with the dot-dot Green 
function. Eventually, the explicit form of the Green's 
functions can be found easily using the EOM technique. 
This way we close the set of mean field equations, which 
are self-consistently solved for each set of parameters (the 
dot levels e^, the tunneling amplitudes V^^i, the flux 0, 
the band width D, and the applied dc bias Vdc)- 

At zero bias we can derive analytical expressions of 
Tk (0-dependent) within the SBMFT. E.g., for 7r-flux 

we get t|"^^^ = i:'exp(-7r|e<j|/2r) (F = t:pq\V\^ is the 
hybridization width). As expected, it is in agreement 
with scaling theory, see Eq. (|2.22fl . 



2. Case Ui2 oo 

For Ui2 ^ cx) only one dot can be charged at a given 
time. In this case we introduce one boson field and one 
constraint that preserves the condition (ni -1-712) = 1. 
The rest of the calculation follows the lines exposed 
above. We find for the Kondo temperature at = tt 
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r: 



^ - (D/^/2)exp(-7r|e<J|/4^) [cf. Eq. ^M]- 



SU(4) 



3. Transport properties 

In the reminder, we describe how to calculate the cur- 
rent through the double QD system within the SBMFT. 
The simplicity of our approach allows us to write the 
current using the Landauer-Biittiker formula: 



, 2e f de ^, 



{!L{e)-fn(e)\ 



(2.41) 



where .'^ie., Vdc) is the transmission probability which de- 
pends on renormalized parameters. Following Meir and 
Wingreen^Si, the transmission through this system can be 
obtained using 



(2.42) 



Here G"'^''^ is the matrix of the advanced (re- 
tarded) Green's function for the dot electrons; i.e., 

G\[rA^) = T«0(±i)({rfz,<.(t),4,.(O)}>- is the ma- 
trix of the renormalized hybridization parameters; i.e., 
ff;ij = TTpiWi^iW*g.jb,b* for Ui2 = and fe-i^j = 

TTpeWe^iWljl]'^ for Ui2 = oo. 

For identical dots and symmetric junctions the trans- 
mission probability is given by 



^(e) = 



id)^ cos^ 



■ 2 (J> 

sm f 



+ {e- SdW^ 



(2.43) 

regardless whether U12 = or U12 — 00. Of course, the 
renormalized coupling F in the above equation should be 
obtained according to the different set of mean-field equa- 
tions, depending on whether U12 = or U12 = 00. We no- 
tice that Eq. (|2.43() was previously obtained in Refsi^Si^ 
for the nonintcracting case. 

The expression for the nonlinear conductance is 
straightforward from the current expression I: — 
dl/dVdc- In the same way, the linear conductance % is 
determined upon insertion of the total transmission eval- 
uated at the Fermi energy into the well known formula 



2e2 



(2.44) 



Following the standard NRG procedure; 



.49.50.51 



we 



evaluate the various physical quantities from the recur- 
sion relation {N > 0) 



AT+l 



N 



^^+1 E E (/Ia^.^/a'.a^+i.- + ^■^) ' (2-45) 



with the initial Hamiltonian given by 



2^ 



1 

71 



^13+ E E^A«Kw + /i-' 



fi—e.o a 



(2.46) 



Here the fermion operators fp,^N,(T have been introduced 
as a result of the logarithmic discretization and the ac- 
companying canonical transformation, A is the logarith- 
mic discretization parameter (we choose A = 2), 



1 - A 



-N 



A-(2N-i)][l- A-(2^'+i)] 



(2.47) 



and 



spectively, are given by 



(2.48) 

The coupling constants Ve and Vo, re- 



K = 4C\/— cos(0/4). 




(2.49) 
(2.50) 



The Hamiltonians J^n in Eq. H2.45|l have been rescaled 
for numerical accuracy. The original Hamiltonian is re- 
covered by 



(2.51) 



with = CA(^-i)/2 . 

In the following, we study the local Green's func- 
tions (with which the linear conductance is calculated) 
and the dynamic spin susceptibility. To improve ac- 
curacy at higher energies, we adopt the density-matrix 
NRG method (DM-NRG)^. In this method, first usual 
NRG iterations are performed down to the energy scale 
lon = -DA~^/^ <^ Tk- From the excitation spectrum at 
this scale, the density matrix is constructed: 



P ■ 



E< 



In 



(2.52) 



C. Numerical Renormalization Group 

SBMFT does not take fully into account real charge 
fluctuation effects. In order to confirm our previous re- 
sults and make quantitative predictions we also use the 
NRG proceduro^'^iSQiSiiS^. 



where |to)^ is the the eigenstate of J^n with energy i?„ 
Then, the NRG iterations are performed again, but now 
at each iteration N' , calculating the Green's function by 



G^a;M'a'(i) = N ' [d A^) ^ d^, J , (2.53) 
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10 12 

V /XSU(4) 
dc K 



FIG. 2: SBMFT results: Left panel {Uu 0): (a) Linear 
conductance (5fo) versus flux for different level positions. 
When the Kondo state is formed (for Ed = —3.5) the 
are delta-like peaks of height 1 centered at even multiples of 
(/)/7r. (c) Curves for 5f versus voltage bias for Ed = —3.5. 
Here we change the flux from (full line) to tt (dot-dot- 
dashed line). Right panel {U12 00): (b) linear and (d) 
differential conductance. Energies are measured in units of 
r = Trpo\V\'^ = D/60. 



where 



PN' = TrAr>7V'/0 



(2.54) 



is the reduced density matrix for the cluster of size N'. 
The Green's function in Eq. (|2.53|) is vahd at the fre- 
quency scale Lo ~ ll>n'- The spin susceptibility is calcu- 
lated in the same manner: 

1 f°° 1 
X{uj) = Im / dt e+'^'-TTpN,{s4t), s,} , (2.55) 



where 



jE[<T'^'^-T-<irfM]- (2-56) 



III. NUMERICAL RESULTS 

We now present our results for the electronic trans- 
port in both limits of the Coulomb interaction, U12 — > 
and U12 00. In the numerical calculations, the 
model parameters are taken as follows: symmetric cou- 
plings (ri(2)L = ri(2)_R = r/2) and equal level positions 
(ei = £2 = Ed)- Throughout this paper, all the parame- 
ters are given in units of the bare coupling F. The energy 
cutoff is set as D = 60r. 



Case U12 



In the left panel of Fig. [21 we present our results when 
C/12 — > obtained with slave-boson mean-field theory. 
First, we focus on the pure Kondo regime when = —3.5 
and discuss both the linear conductance and the non- 
linear conductance (given by % = ^^(Vdc = 0) and 
= dl/dVdc, respectively). The linear conductance 
[solid line in Fig. I2a)] shows narrow peaks due to con- 
structive interference around « (mod 2tt) whereas 
transport is suppressed elsewhere. This is caused by the 
fact that the DOS of each dot has a resonance exactly at 
Ep. In the language of slave bosons this means £1^2 = 
and a SU(2) Kondo state is well formed. Therefore, these 
narrow peaks in '^q correspond to paths through the AB 
geometry with multiple windings around the enclosed 
flux. The width of each peak is given roughly by « Tk- 
Away from the constructive interference condition, the 
transmission at the Fermi energy quickly vanishes. This 
unusual behavior is clarified with our calculations of the 
differential conductance. In Fig.|2Ic) we show the nonlin- 
ear conductance dl/dVdc as a function of the bias voltage 
Vdc for different values of the flux and £;; = — 3.5. In the 
absence of flux (or for even multiples of it) the nonlinear 
conductance shows the usual zero-bias anomaly (ZBA), a 
narrow peak at Vdc = that reaches the unitary limit due 
to the constructive interference in the resonant condition. 
Increasing (j> does not affect the Kondo resonance much 
so that the transmission probability can be written 
as a combination of a Breit-Wigner resonance for id = 
plus a Fano antiresonancei^S A dip at zero bias is then 
obtained [see Fig. |2Ic)]. The width of this dip is 

Tr-(0)[1 — cos(20)]. It has an oscillatory dependence on 
the applied flux. This result is in good agreement with 
the NRG calculations as shown in Fig. |3fa). Here, we 
plot ^ as a function of energy. It is worthwhile to note 
that ^ amounts to at low bias. 

For increasing Sd one enters the mixed-valence regime 
[see Fig. |2Ia)]. Although the results should be taken 
in a qualitative way, we find that the renormalized lev- 
els for Sd — —2, —1.75, —1.5 are no longer at Ep except 
when (/) « (mod 27r). The transmission coefficient (and 
thereby the conductance) is extremely sensitive to devi- 
ations of id out of Ep. When the bare level position is 
shifted toward the Fermi energy the renormalized levels 
for Ed — —2, —1.75 as a function of are not at Ep except 
when w (mod 27r) whereas for Sd = —1.5 they never 
reach Ep. In these cases due to the lack of a resonant 
condition at each dot, multiple windings are less likely to 
occur and the conductance starts to resemble a cosine-like 
function generated by a combination of lower harmonics 
[see Fig. |2a), case Sd = —1.5]. For Sd = —2,-1.75 we 
still observe the sharp resonance at $^(0 « 0) (mod 27r) 
due to a quasiresonant condition when (f> ^ (mod 27r). 
Finally, for Sd = —1.5 the linear conductance has a trivial 
cosine dependence. 
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B. Case U12 00 

Next, we elaborate on the numerical results for the 
limit of a strong interdot Coulomb interaction [/12 ^ 00 
(right panel of Fig. [^J. The results show that in this 
situation not only the spin fluctuates but also the pseudo- 
spin since just two charge states are allowed in the double 
QD system: {ni = \,n2 = 0} and {ni = 0,n2 = 1}- The 
fluctuations in both sectors (spin and pseudo-spin) leads 
to the exotic SU(4) state close to (/> = tt (mod 2tt). 

We begin with the linear regime. Figure |2Ib) sum- 
marizes our results for as a function of the applied 
flux. We concentrate on the pure Kondo regime and set 
£d = —7 well below Ep- Unlike the case of weak interdot 
Coulomb interaction [see Fig. [^Ja)], the linear conduc- 
tance shows broad peaks at positions w (mod 27r). In 
addition, the linear conductance only vanishes when the 
condition of destructive interference takes place. Let us 
investigate in in some detail the two limit cases « 
and « TT (mod 2n). In our RG analysis, we find for 
« (mod 2ti) that the ground state corresponds to the 
usual spin SU(2) Kondo effect with a greatly enhanced 
Kondo scale. Accordingly, the corresponding renormal- 
ized level lies at = leading to a shift of the scattering 
phase 5 = 7r/2. On the contrary, for cj) — tt (mod 27r) we 
find that the ground state is the highly symmetric SU(4) 

Kondo state with a renormalized level at id ~ Tfy^'^\ 
which implies (5 = 7r/4 to fulfill the Friedel-Langreth sum 
rulei^ Quite generally, in a SU(N) problem the phase shift 
becomes S = tt/c/K in the limit of large ^ and the Kondo 
resonance shifts up to id ~ ttF/o/K.^ In the intermediate 
regime, when < < tt, the renormalized level takes 
on a positive value id < T^'"'^\ As a consequence, away 
from « (mod 2tt) the resonant condition is not sat- 
isfied (the renormalized level id is not longer at Ep). 
In this situation electronic paths with multiple windings 
do not occur and the linear conductance consists of a 
cosine-like function [Fig. Efb)] . At finite Vdc the nonlin- 
ear conductance displays a ZBA which is quenched as 
(j) decreases [see Fig. Eld)], unlike the dip found in the 
U12 = case. Eventually, for = tt (mod 2tt) there is no 
transport due to completely destructive interference. 

We can compare our results shown in Fig. EJd) with 
those obtained from NRG plotted in Fig. |3Ib). Here, 
one can see that =5^ decreases as </> increases, which is 
consistent with the results of the SBMFT. Nevertheless, 
SBMFT overestimates the decreasing rate of the ZBA. 
The NRG results show that while the peak does not 
change appreciably for (/> < 0c, it decreases very rapidly 
for d) > (hr- 



C. Crossover 

The value of 0c is the last ingredient we have to ex- 
plain. 0c marks the crossover between SU(2) Kondo 
physics to the highly symmetric SU(4) Kondo state. For- 
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FIG. 3: NRG results: Top panel: Transmission probability 
versus flux for (a) Ed = — 7r, Ui — U2 = 5D, U12 ~ 0, and 
(b) £d = -14r, Ui = U2 = 5D, and U12 = 5D. We set 
r — D/60. (Notice that we do not recover the unitary limit 
of 5^ for ~ (mod 2-k) because of the systematic errors 
introduced in the NRG procedure). Bottom panel {U12 = 
5D): (c) Spin susceptibility (in an arbitrary unit) in the limit 
of strong interdot interaction, (d) The peak position of the 
susceptibility as a function of the flux 0. 



tunately, 0c can be extracted from the peak position 
of the spin susceptibility xi^)i which yields a reason- 
able estimate of the Kondo temperature. Figure |3fc) 
shows the evolution of x when increases. Remarkably, 
when the flux enhances, at some point the position of 
the peak moves toward higher frequencies. The peak 
position as a function of is plotted in Fig. |2td) . We 
observe that Tk (0) is almost constant when goes from 
zero to 0c ~ 0.757r. This fact allows us to establish a 
critcrium for the crossover between the SU(2) and SU(4) 
Kondo states in the double QD system. 



IV. CONCLUSION 

We have analyzed the transport properties (in and out 
of equilibrium) of a prototypical mesoscopic double-slit 
interferometer when interactions play a dominant role. 
We have shown that crucial differences arise in the lim- 
its of negligible and large the interdot Coulomb interac- 
tions. In the former case, only spin fluctuations matter 
and each dot develops a Kondo resonance at the Fermi 
level independently of the applied magnetic flux. Due to 
the interference between these two Kondo resonances the 
linear conductance versus the flux shows a series of nar- 
row peaks at = 27r (mod 2tt) of unitary height (in units 
of 2e^//i). Furthermore, we have found that any devia- 
tion from the Kondo regime (close to the mixed-valence 
regime) lead to dramatic changes in the conductance as 
a function of the flux. Interestingly, the nonlinear con- 
ductance show the formation of a dip when 7^ 27r (mod 
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27r). A complete suppression of the electronic transport 
occurs when the destructive interference condition takes 
place (j)Tr (mod 2tt). 

Charge and spin become entangled when the inter- 
dot Coulomb interaction is very large. Here, the dif- 
ferential conductance has a zero bias anomaly quenched 
with increasing flux. The Kondo state changes its sym- 
metry, from SU(2) to SU(4), as (j) approaches tt (mod 
2tt). Since the crossover is not too close to cf) — tt 
(mod 27r) the SU(4) state remains robust to be detected 
experimentally. Our geometry requires symmetric cou- 
plings to the leads but not inevitably equal (i.e., we need 
^iL + T2R = riK-Fr2L).SThe charging energies Ui, U2, 
and U12 should be of the same order (a few meV). Finally, 



the external flux should correspond to a low magnetic 
field to avoid spin Zeeman splittings in the dot, around 
10 niTi2£ All these constraints are experimentally acces- 
sible with present techniques 
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